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1. Introduction 



The scope of this article is two-folded: the first part is dedicated to the analysis of 
the propagation uncertainties of cosmic rays in the Milky Way; in the second part we 
studied the detection possibilities of an eventual neutralino induced component in the 
antiproton spectra in the framework of minimal supergravity (mSUGRA). 

In order to disentangle an exotic contribution from the standard one in the 
antiproton spectra first we have to test the accuracy of the standard calculation. 

In Section 2 we give a brief description of standard mechanisms for the 
propagation of cosmic rays and we present the results about the uncertainties 
in the standard propagation model due to the unknown geometrical and 
hydromagnetodynamical parameters of the Galaxy, uncertainties of the measurements 
and parametrizations of nuclear cross sections and uncertainty of the helium to 
hydrogen ratio in the Galaxy. 

In Section 3 we describe creation and propagation of a neutralino induced 
component. 

In Section 4 we study the possibility to detect secondary components in the 
positron and antiproton fluxes by the upcoming PAMELA satellite experiment. 

In Section 5 we present the results about the possibility to disentangle an eventual 
signal component in the antiproton spectra. We did this analysis in a clumpy halo 
scenario. Finally we found the minimal and the maximal values of the dumpiness 
factors needed to disentangle the signal from the background without violating present 
antiproton data. We also give some examples of total fluxes in comparison with the 
experimental data. 

Section 6 contains the results and the conclusions. 



2. Propagation of Cosmic Rays in the Milky Way and Its Uncertainties 

The most complete equation for the propagation of cosmic rays that includes all the 
known physical processes is 
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where ip(r,p,t) is the total phase space density. We will shortly review here the 
main features of the physical processes described by this equation implemented in the 
Galprop code 

The second term describes isotropic diffusion, defined by the coefficient that 
depends on the rigidity (momentum per unit of charge, p = p/Z) 

D xx = f3D ( P /po) S , (2) 

inspired by the Kolmogorov spectrum (8 = 1/3) of the weak magnetohydrodynamic 
turbulence. In [3] was first shown that the Kolmogorov spectrum best reproduces the 
sharp peak in B/C data. In some models we used a break in the index S at some 
reference rigidity p : Si = below p while 5 2 > above p . 

The convection velocity field V c , that corresponds to the Galactic wind, has 
a cylindrical symmetry. Its z-component is the only one different from zero. It 
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increases linearly with the distance z from the Galactic plane. This is in agreement 
with magnetohydrodynamical models 0]. In the Galactic plane there should be no 
discontinuity in the convection velocity field and so we considered only V c (z = 0) = 0. 

Reacceleration is determined by the diffusion coefficient in the momentum space 
D pp . D pp is a function of the corresponding configuration space diffusion coefficient 
D xx and of the Alfven velocity Va in the framework of quasi-linear MHD theory 

nana ei 



4jW 

36(4 - 5 2 )(4 ~ 6)w ' 



Dp P {D XXl V A ) = oUA P „J A , (3) 



where w characterizes the level of turbulence, and it is equal to the ratio of MHD 
wave energy density to magnetic field energy density. It is assumed w = 1, but the 
only relevant quantity is V\jw. 

The unknown values of parameters such as the Alfven velocity, the convection 
velocity gradient in Milky Way and the height of the galactic halo can be constrained 
by the B/C data. With the sets of the constrained parameters one can find all 
the possible spectra for the others cosmic rays. This procedure was already used 
in for another propagation code. 

Injected spectra of all primary nuclei are power laws 

dq(p)/dpocp" 1 , (4) 

where 7 is the injection index. In principle its value can vary with species. It can 
be shown that the power law approximation as well as a small break in the injection 
indexes 7 is allowed in the framework of diffusive shock acceleration models |lllll2llT3*| . 

Source term q(r,p) for secondaries contains cross sections for their production 
from progenitors on H and He targets 

q{r,p) = fJcipp(r,p)[o-^(p)n H (r) + a p ^ e {p)n He (r)], (5) 

where a p ^ (p) and a p ^ e (p) are the production cross sections for the secondary from the 
progenitor on H and He targets, ip v is the progenitor density, and n#, are the 
interstellar hydrogen and helium number densities. 

The last two terms in equation Q are loss terms with characteristic times 77 and 
t> related respectively to the fragmentation and radioactive decay. The heliospheric 
modulation of the local interstellar spectra in the vicinity of the Earth and in the 
heliosphere hole has to be taken into account in order to obtain the realistic cosmic 
rays spectra in locations where they are/will be measured (balloon-born or satellite- 
borne experiments). We made use of a widely used and tested model in which the 
transport equation is solved in the force field approximation |141 [To] . That equation 
describes diffusion processes in the heliosphere and includes effects of heliospheric 
magnetic field and solar wind. In this case, solar modulation is a function of just 
a single parameter that describes the strength of the modulation. All the dynamical 
processes are simulated simply changing the interstellar spectra during the propagation 
inside the heliosphere: 
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Figure 1. Propagation parameters uncertainty for B/C ratio: for the DR model 
it is given with solid lines around the best fit (dashed line), while for the DC 
model it is given with dotted lines around the best fit (dashed line). For DRB 
model we give the best fit (dashed line). For the complete list of the experimental 
data see |f 6| . 



where E and p are energies and momenta of the interstellar and of top of the 
atmosphere fluxes. The solar modulation is entirely determined by the <j) parameter. 

In Galprop the model of the Galaxy is three dimensional with cylindrical 
symmetry; the coordinates are (R,z,p), where R is Galactocentric radius, z is the 
distance from the Galactic plane, and p is the total particle momentum. The distance 
from the Sun to the Galactic centre is taken to be 8.5 Kpc. The propagation region 
is bounded fixing R max = R = 30 Kpc and z max = z beyond which free escape 
is assumed. The distribution of cosmic rays sources is chosen to reproduce (after 
propagation) the cosmic rays distribution determined by the analysis of EGRET 
gamma-ray data done in 13 . The code first computes propagation of primaries, 
giving the primary distribution as a function of (R,z,p). Then the secondary source 
function is obtained from the gas density and cross sections. Finally, the secondary 
propagation is computed. 

Secondary to primary CR ratios are the most sensitive quantities on variation of 
the propagation parameters. This can be verified numerically. Primary to primary 
and secondary to secondary ratios are not very sensitive to changes in the propagation 
parameters because they have similar propagation mechanisms. The most accurately 
measured secondary to primary ratio is boron to carbon ratio (B/C, see JE]). Boron 
is secondary while, one of its progenitors, carbon is primary. The B/C data are used 
also because they have relatively well known cross sections. To estimate the quality 
of the data fit we used the standard \ 2 test 

Y 2 = 1 y \ ($b/c „$s/c )2 (8) 

n (°~n ) 

where a B / c are statistical errors for N = 46 experimental points and $>elp are 
measured and &fJ C are predicted values of the ratio. 
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Table 1. Allowed values for diffusion and reacceleration model propagation 
parameters. 



par./val. 


z [Kpc] 


D [cm 2 s x ] 


S 


7 


va [Kms x ] 


minimal 


3.0 


5.2 10 2H 


0.25 


2.35 


22 


best fit 


4.0 


5.8 10 2H 


0.29 


2.47 


26 


maximal 


5.0 


6.7 10 28 


0.36 


2.52 


35 



We analyzed only the two extreme cases of the propagation models: either without 
convection or without reacceleration. We did not succeed in obtaining satisfactory 
models with all the physical processes switched on. A possible explanation is that 
one of the processes is really subdominant in the Galaxy in such a way that the 
data require either reacceleration or convection to be well fitted. So in the first case 
we included diffusive and reacceleration effects (see equation QJ and we considered 
two sub-cases: (DR) without a break in the index of primary injection spectra and 
(DRB) with a break at the rigidity p 1 . The second case (DC) contains diffusion and 
convection terms from the same equation |T]l. In this model there are two breaks: 
one in the index of the primary injection spectra and the other in the spectra of the 
diffusion coefficient D xx . The first break means that at some rigidity p 7 the index 
7 suffers a discontinuity (see equation and ref. f° r details), while the second 
one means that at some rigidity po the diffusion index 6 suffers a discontinuity (see 
equation (0) and discussion in ref. 0). 

For the DR model we chose to vary the following parameters: the height of the 
Galactic halo z, the constant in the diffusion coefficient Do (from equation the 
index of the diffusion coefficient 6 (from the same equation), the primary spectra 
injection index 7 for all the energies (from equation ijljl) and the Alfven velocity va 
that determines the strength of reacceleration. In order to find the allowed range 
of these parameters we have required a reduced x 2 l ess than 2 for the fit of the 
B/C experimental data [T^]. In figure ^ are presented the enveloping curves of all 
the good fits with solid lines around the best fit line for the same model, that is 
represented with dashed line. We took the experimental data with relatively small 
solar modulation parameter cf> between 325 MV and 600 MV, where the force field 
approximation is better justified than for higher modulation parameters. The allowed 
ranges of propagation parameters of this model are given in the table Q Using this 
result we found the enveloping curves of all the positron and antiproton spectra. In 
this way we determined the upper and the lower bounds of the uncertainty bands for 
positron and antiproton spectra due to the propagation parameters uncertainty (see 
figure [21 for the positron spectra uncertainty and figure [3] for the antiproton spectra 
uncertainty). The relative uncertainty depends on the energy range of the spectra. 
For positrons the relative uncertainty varies from 30% under 1 GeV to 15% around 
10 GeV, increasing again after 10 GeV (figure 0) while for antiprotons it varies from 
about 10% up to about 15% (figure 0). 

For a better visualization some figures that contain positron and antiproton data 
contain just the part of the experimental data chosen to cover in a convenient way all 
the energy range reached in experiments. The data are taken from references |19II17| . 

In the case of the DC model we chose to vary the following parameters: Dq, 
the diffusion index 5% below the reference rigidity po = 4 GV and 62 above it (see 
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Figure 2. Upper and lower bounds of positron spectra due to the uncertainties 
of the propagation parameters are represented with solid lines for the DC model, 
with dashed lines for the DR model and with dotted lines for the DRB model. 
The DR and DRB model uncertainties are very similar, but there is a slight 
improvement of the fit in the low energy part of the spectra in the case of DRB 
model. On the other hand, around the maximum the DRB model overestimates 
the data slightly more than the DR model. Experimental data are taken from 1171 . 



DC model 
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Figure 3. Upper and lower bounds of the antiproton spectra due to the 
uncertainties of the propagation parameters for the DC model (solid lines) and 
DRB model (dotted lines). The propagation parameters uncertainties for the DR 
model are almost the same as those of the DRB model (see the discussion later in 
the text), so we are not presenting them here to avoid the confusion. Experimental 
data are taken from I19| . 
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Table 2. Allowed values for the propagation parameters for diffusion convection 
model. 



par./val. 


z [Kpc] 


D [s£] 


s 2 


dVc r Km i 
dz L skpc -1 


7i 


72 


minimal 


3.0 


2.3 10 28 


0.48 


5.0 


2.42 


2.14 


best fit 


4.0 


2.5 10 2H 


0.55 


6.0 


2.48 


2.20 


maximal 


5.0 


2.7 10^ 


0.62 


7.0 


2.50 


2.22 



Table 3. Allowed values for the propagation parameters for diffusion 
reacceleration model with break in the injection spectra of primary nuclei. 



par./val. 


z [Kpc] 




5 


7i 


72 




minimal 


3.5 


5.9 10^ 


0.28 


1.88 


2.36 


25 


best fit 


4.0 


6.1 10 28 


0.34 


1.92 


2.42 


32 


maximal 


4.5 


6.3 10^ 8 


0.36 


2.02 


2.50 


33 



equation J2J), the halo size z, the convection velocity V c (see equation JU) and the 
injection index for primary nuclei 71 below the reference rigidity = 20 GV and 72 
above it (see equation Q)). Enveloping curves of the B/C fits for a reduced x 2 values 
less than 2.8 are given in figure^ Positive variations around Si — gave unsatisfactory 
fit. In order to take the smallest possible break of this index we decided not to take 
negative 5± values. Allowed values for the propagation parameters can be found in 
table |21 Again, we used them to derive the uncertainties of positron and antiproton 
spectra. Relative uncertainty for positrons vary between 20% above the maximum 
and 30% below it (figure [2J while for antiprotons is about 20% around 20 MeV, 17% 
around the maximum and 25% around 20 GeV (figure EJ. 

We also found the spectra that correspond to the parameters of the best fit of 
B/C data for subFe/Fe ratio (see figure 0}, protons, helium and electrons as well as 
the corresponding propagation parameters uncertainties. The solar activity (around 
a minimum in 1998 and around a maximum in 2002) is the cause of the discrepancies 
in the data below 10 GeV for protons. For DC model the fits are good, while DR 
overestimates protons (figured), helium (figure EJ) and electrons (figure EJ). 

In order to improve those fits, we considered the DR model with a break in the 
injection index for the primary nuclei spectra with a rigidity of 10 GV [21 E3- We 
determined the allowed values of the propagation parameters (table |SJ demanding the 
same reduced x 2= 2 as for the DR model (see figure 0). The positron and antiproton 
uncertainties are presented in figures [21 and Even if positrons at low energies 
and protons and helium in all the energy range gives a better fit (see figure and 
figure EJ), they remain overestimated. For the computation of B/C ratio, Galprop 
uses only one principal progenitor and compute weighted cross sections. Introducing 
the break in the index of the primary injection spectra in the DR model give worst 
electron data fit (see figure EJ than in the case without the break. The antiproton 
spectra remain unchanged, still significantly and systematically underestimated in all 
the energy range. 

We also calculated the variation of the antiproton spectra as a function of the 
most important antiproton production cross sections. Antiprotons are created in 
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Figure 4. Ratio (Sc+Ti+V) /Fe that corresponds to the propagation parameters 
that give the best fits of B/C data for the DC model given with dashed line. It 
is inside the corresponding uncertainty band also given with dashed lines. The 
ratio for the DR model is given with dotted line and it is inside the uncertainty 
given with solid lines, while for DRB model is given with larger-step dashed line 
without the uncertainty band around. Experimental data are taken from 121)1 . 




Figure 5. Top of the atmosphere spectra of electrons that correspond to the 
parameters of the best B/C fit are the lower curves: for the DC model they are 
given with solid line, for the DR model with dashed line and for DRB model with 
dotted line. The local interstellar spectra are the upper curves; the three models 
are represented with the same types of lines. Experimental data are taken from 

El 
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Figure 6. Upper and lower bounds of helium spectra due to the uncertainties of 
the propagation parameters for the DC model are represented with dashed lines. 
Spectra that correspond to the parameters of the best B/C fit are given for the 
DC model with solid line, for the DR model with dashed line and for the DRB 
model with dotted line. Experimental data are taken from 1181 . 

the interactions of the primary cosmic rays (protons and other nuclei) with the 
interstellar gas. Dominant processes are interactions of high energy primary protons 
with hydrogen, p + p^p + p + p + p. Parametrization of this cross section is given in 
|21| . Other cross sections, those of the primary protons with other nuclei, are studied 
in reference |22j . From these, the most important are those that involve helium, and 
they contribute less than 20% of the total production of all the antiprotons. All the 
heavier nuclei together give just a few percents of the total production. 

Uncertainties of cross sections influence the antiproton spectra uncertainties: 
simultaneous settings of all the production cross sections to the maximum/minimum 
rise/lower the upper/lower propagation parameters uncertainty bounds (already found 
before). Errors obtained in this way give contributions to the total uncertainties: they 
vary from 20% up to 25% in the case of a DR model and, very similarly, from 20% up to 
24% for a DC model (depending on the energy range of the spectra). Non-production 
cross sections, the so called tertiary component, correspond to inelastically scattered 
secondaries p+X — > p+X. Those processes bring down the energies of the antiprotons 
of relatively high energies, flattening like that the spectra. But, even if the uncertainty 
related to those cross sections is relatively big, this does not give relevant change of 
the antiproton spectra because tertiary contribution is very small. In fact, it has been 
implemented in propagation codes just recently. 

The uncertainty in the measurements of helium to hydrogen ratio brings another 
component to the total uncertainty of the cosmic rays spectra. By changing the Hc/H 
ratio in a reasonable range from 0.08 to 0.11 (see and [23) we obtained a relatively 
small contribution. For both positron and antiproton flux uncertainty and for all of 
the models we considered it varies from 3% to 7%, depending on the energy. 
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Figure 7. Upper and lower bounds of proton spectra due to the uncertainties of 
the propagation parameters for the DC model are represented with dashed lines. 
Spectra that correspond to the parameters of the best B/C fit are given for the 
DC model with solid line, for the DR model with dashed line and for the DRB 
model with dotted line. Experimental data are taken from 1181 . 



Total uncertainties of positrons and antiprotons are presented in figure and in 
figure [5] respectively. They vary from 35% up to 55% for antiprotons and from 20% 
up to 40% for positrons for both the models in the current experimental data energy 
range. 

In plot 1101 we presented two different uncertainties of the positron spectra due 
to the solar modulation for the same DRB model (that fits better the data than the 
DR model) . The first uncertainty is obtained solar-modulating the lower bound of the 
propagation uncertainty band using the 10% bigger <f> than the measured one and the 
upper bound with 10% lower modulation parameter. This is a conservative error for 
the solar modulation parameter. The spectra are still completely above the data. As 
we used the force field approximation for the solar modulation, we tried also to change 
the modulation parameter for ± 50%. The resulting curves are represented on the plot 
and are still completely above the data. As we deal with a medium solar modulation 
(4> 600MV) for positron data (we are not near the solar minimum) the sign-charge 
drift effect due to the solar magnetic field polarity should be not very strong. 

3. Component of the Antiproton Spectra Induced by Neutralino 
Annihilations 

In this section we take into account the possibility of a neutralino induced component 
in the p flux. Our analysis is performed in the well known mSUGRA framework 25 
with the usual gaugino mass universality at the grand unification scale Mqjjt- 

In the general framework of the minimal supersymmetric extension of the 
Standard Model (MSSM), the lightest neutralino is the lightest mass eigenstate 
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Figure 8. Total uncertainties of positron fluxes and spectra that correspond 
to the parameters of the best B/C fit for the DC model (solid lines around the 
best fit curve, also solid), DR model (dashed lines around the best fit curve, also 
dashed) and the DRB model (dotted lines around the best fit curve, also dotted). 
Experimental data are taken from |17| 

obtained from the superposition of four interaction eigenstates, the supersymmetric 
partners of the neutral gauge bosons (the bino and the wino) and Higgs bosons (two 
Higgsinos). Its mass, composition and couplings with Standard Model particles and 
other superpartners are functions of several free parameters one needs to introduce 
to define such supersymmetric extension. In the mSUGRA model, universality at 
the grand unification scale is imposed. With this assumption the number of free 
parameters is limited to five 

m i/2i m 0i sign(fi), Ag and tan/3, 

where mo is the common scalar mass, is the common gaugino mass and Aq 

is the proportionality factor between the supersymmetry breaking trilinear couplings 
and the Yukawa couplings, tan/3 denotes the ratio of the VEVs of the two neutral 
components of the SU(2) Higgs doublet, while the Higgs mixing /j, is determined (up 
to a sign) by imposing the Electro- Weak Symmetry Breaking (EWSB) conditions at 
the weak scale. In this context the MSSM can be regarded as an effective low energy 
theory. The parameters at the weak energy scale are determined by the evolution 
of those at the unification scale, according to the renormalization group equations 
(RGEs) HI]. 

For this purpose, we have made use of the ISASUGRA RGE package in 
the IS A JET 7.64 software [2Z| - After fixing the five mSUGRA parameters at 
the unification scale, we extract from the ISASUGRA output the weak-scale 
supersymmetric mass spectrum and the relative mixings. Cases in which the lightest 
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Figure 9. Total uncertainties of antiproton fluxes and spectra that correspond 
to the parameters of the best B/C fit for DC (solid lines around the best fit curve, 
also solid) and DRB model (dotted lines around the best fit curve, also dotted). 
Experimental data are taken from 1191 



1000 




Figure 10. Top of the atmosphere (lower solid curve) and local interstellar 
(upper solid curve) spectra of positrons that correspond to the parameters of the 
best B/C fit for DRB model. Total uncertainties of positron fluxes are presented 
with dashed lines around the best fit. Two enlargements of the propagation 
uncertainties due to two different solar modulation uncertainties are presented 
with dotted and dashed lines (see the text for details). Experimental data are 
taken from I17| . 
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neutralino is not the lightest supersymmetric particle or there is no radiative EWSB 
are disregarded. 

The ISASUGRA output is then used as an input in the DarkSUSY package j.''>(i|. 
The latter is exploited to: 

• reject models which violate limits recommended by the Particle Data Group 2002 
(PDG) EH; 

• compute the neutralino relic abundance, with full numerical solution of the 
density evolution equation including resonances, threshold effects and all possible 
coannihilation processes 

• compute the neutralino annihilation rate at zero temperature in all kinematically 
allowed tree-level final states (including fermions, gauge bosons and Higgs 
bosons); 

• DarkSUSYestimates the induced antiproton yield by linking to the results of the 
simulations performed with the Lund Monte Carlo program Pythia j.'iOj . 

This setup as well as some other similar scenarios were already considered in the 
context of dark matter detection and of an improvement of the cosmic rays data fits 
(a list of references includes, for example, |311 \'32\). The comparison of our results 
with previous works and other complementary techniques should be transparent. 



3.1. Clumpy Halo Models 

The dependence of the antiproton flux is oc p 2 /ra 2 (see H10|) . In the case of a small 
clump scenario |33| for the dark matter halo in our Galaxy there is a local density 
enhancement without increasing the total halo mass and then an increase of the 
antiproton flux. 

By hypothesis the clump is a spherical symmetric compact object with mass 
M c i and some density profile p c i (r c i). We denote with / the dark matter fraction 
concentrated in clumps and we introduce the dimensionless parameter d 

d= 1 /d 3 r cl [p cl (r cl )] 2 
Po J d?r c \p c \ (f c i) 

that gives the overdensity due to a clump with respect to the local halo density 
po = p(ro), where ro is our distance from the Galactic Center (GC). In a smooth 
halo scenario the total neutralino induced p flux calculated for r = r is given by PI] 

$,(r , r) ee (<w) £ 2 C prop (T) . (10) 

f \ m xJ 

where T is the p kinetic energy, a ann v is the total annihilation cross section times 
the relative velocity, m x is the neutralino mass, B' and dN* /dT, respectively, the 
branching ratio and the number of p produced in each annihilation channel / per unit 
energy and C pr0 p(7 1 ) is a function entirely determined by the propagation model. In 
the presence of many small clumps the p flux is given by 

*f naw (r ,T)=fd-$ p {r ,T) (11) 

For the smooth profile we assumed a Navarro, Frenck and White profile (NFW) |35| . 
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3. 2. Propagation of the Neutralino Induced Component 

The primary contribution to the antiproton flux is computed using the public code 
DarkSUSY[37]. We modified the antiproton propagation in order to be consistent with 
the DC propagation model as implemented in Galprop code. We assumed diffusion 
coefficient spectra used in Galprop code with our best fit values for the diffusion 
constants Do and 6. In DarkSUSY, the convection velocity field is constant in the upper 
and the lower Galactic hemispheres (with opposite signs, and so it suffers unnatural 
discontinuity in the Galactic plane), while Galprop uses magnetohydrodynamically 
induced model, in which one component of velocity field along the Galactic latitude 
(the only one that is different from zero) increases linearly with the Galactic 
latitude 0] . We assumed an averaged convection velocity calculated from the Galactic 
plane up to the Galactic halo height z. 

4. Detection of the Secondary Components in the Positron and the 
Antiproton Fluxes by PAMELA 

In this section we calculate the statistical errors for the PAMELA |H81 0Hj experiment 
for the positron and the antiproton background spectra calculated with Galprop L 
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Figure 11. Geometric factor of PAMELA 

The calculation is done for a three years mission assuming the geometric factor 
given in figure 1111 taken from the results of the test measurements on triggers with 
the first and the last scintillators of the apparatus @D]. PAMELA'S statistical errors 
for positrons and antiprotons in the case of the DC model are given in figure 1121 and 
in figure We present in figures El an d El PAMELA expectations together with 
the propagation uncertainties for the antiproton proton ratio and positron charge 
fraction, respectively. For completeness we also present the expectations for the DRB 
propagation model for the antiprotons in figure ITU We do not present PAMELA'S 
expectations for antiprotons in the DR model case because the DR model flux that 

| The list of the people and the institutions involved in the collaboration together with the on-line 
status of the project is available at http://wizard.roma2.infn.it/ 



Uncertainties of Cosmic Ray Spectra and Detectability of Antiproton 



15 



10 r 




HEAT 94 
a MASS 89 

■ DC - PAMELA 3 years 

i_j i i i i J i i i 

10" 1 1 10 10 2 

kinetic energy (GeV) 



Figure 12. Experimental data (from 1171 ') confronted with PAMELA'S 
expectations for positrons for DC model background. Total uncertainties of 
positron fluxes with the spectra that correspond to the parameters of the best 
B/C fit in the middle are given for a better comparison: for DC model they are 
represented with dotted lines while for DRB model with dashed lines. 



correspond to the best B/C fit propagation parameters is almost identical with that 
of DRB model. DR and DRB model background PAMELA predictions for positron 
fluxes are not given because those fluxes largely overestimate the experimental data 
at low energies. 



5. The Possibility of Disentanglement of the Neutralino Induced 
Component in the Antiproton Flux With PAMELA 

In this section we present the results we found about the minimal values of the 
dumpiness factors fd needed to disentangle a neutralino induced component in the 
antiproton flux with PAMELA. We computed this factor as a function of the mSUGRA 
parameters, fixing Aq, tan/3 and sign(/i) = +1. In this way the dumpiness factor 
become a function of mo and m^ 2 parameters. Similar analysis were already made 
in the literature (see for example |32)L 

For the discrimination we requested the following conditions: 

(i) The total antiproton flux <j> tot — <f>bkg + (f'susy gives a good fit of the experimental 
data. 

(ii) Difference between </> tot and the DC model <fibkg is detectable by PAMELA. 
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Figure 13. Experimental data (from 1191 ) confronted with PAMELA'S 
expectations for antiprotons for DC model background. Total uncertainty of the 
DC model antiproton background flux with the spectra that correspond to the 
parameters of the best B/C fit in the middle are given for a better comparison; 
they are represented with dotted lines. The reduced x 2 is 



The first condition is satisfied if 

i . . (cb ex P — <J> tot "| 2 

v 2 - 1 V ^ n ) 2 

x f il ~ N _ 1 2^ ( a exp ) 2 fit '° [ 1 

n \ n J 

where x}u o = f° r N ~ 40 experimental points. The second condition is satisfied 
if 

2 1 r (^ 9 -^) 2 > 3 (13) 

Adiscr i / > , Pbka\n — A,discr,0 \- LU / 

M ~ 1 m ( a m ) 2 

where Xdiscr o = *-8> f° r ^ = 2 ^ points and where it^* 9 are the PAMELA statistical 
errors associated to the background flux (those presented in figure IT3)) . 

The reduced x}it f° r ^ nc background flux alone is Xfu = This means that 
the background already gives a good fit of the experimental data. 

5.1. Results 

For each model we found the minimal value of the dumpiness factor fd needed to 
satisfy both conditions. As the dumpiness factor is a function of too and mwj 
parameters we made contour plots calculating equi-clumpincss factors lines. We also 
found the maximally allowed values for fd in order not to violate the condition Ijl2|l . 
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Figure 14. Experimental data (from 1191 ) confronted with PAMELA'S 
expectations for antiprotons for DRB model background. Total uncertainty of 
the DRB model antiproton background flux with the spectra that correspond to 
the parameters of the best B/C fit in the middle are given for a better comparison; 
they are represented with dotted lines. 



First, we found the results for different values of tan/3, presented in figure 1171 and in 
figure 1211 Then we examined in more details some interesting zones in those plots. 
The regions in parameter space that are excluded either by accelerator bounds or 
because electroweak symmetry breaking is not achieved or because the neutralino is 
not the lightest supersymmetric particle are represented with black color. Red (dark 
shaded) color indicates the (mo, mi/2) domains with ilh 2 in the WMAP [JI] region 
0.09 < ilh 2 < 0.13. Green (light shaded) color indicates the parameter space regions 
with values of 0.13 < Vlh 2 < 0.3. Equi-clumpiness factors lines for the minimal fd 
needed to disentangle a supersymmetric signal with PAMELA are given with black 
solid lines. Beside every line it is indicated the value of the dumpiness factor. Under 
the line fd — 1 there are no models that satisfy at the same time the both conditions 
from equations i|12|) and l|13|) . As an example we give in the right panel of figure ITTjl 
a plot with that region colored in blue. In figures El the equi-clumpiness factors lines 
for the maximal allowed (by the present experimental data) fd are represented by the 
lower bound of the translucent light blue (very light shaded) regions. In these figures 
the minimal fd are instead the upper bounds of the same regions. The translucent 
regions denote the parameter space domains that correspond to models that satisfy 
both the conditions (fT2f and lfLl|l , but now keeping fixed values of fd (those indicated 
inside the regions). 

It can be seen that relatively small dumpiness factors (of order 10) are sufficient 
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Figure 15. Experimental data (from I19| and |18l ) confronted with PAMELA'S 
expectations for the antiproton proton ratio for the DC model background. The 
propagation uncertainty band of the antiproton proton ratio and the curve that 
corresponds to the parameters of the best B/C fit in the middle are given for a 
better comparison. 



for PAMELA detection. This is very important, because, even if we consider a DC 
model as the background flux (that alone already gives a good fit of the experimental 
data) it is still possible to disentangle a supersymmetric component in a wide region of 
the parameter space (in comparison with the WMAP allowed zone). The two favourite 
cases are: tan/5 = 55 (see figure an d tan/? = 60 (see figure I2T|) . 

In figure ITH1 we presented one of the two particularly interesting regions, namely 
those inside the cosmologically allowed zone, found for the tan/3 = 50 case, i.e. that 
from the upper panel of figure 1171 The other interesting region, corresponding to 
the so called focus point region 021 m which the neutralino has a significant higgsino 
component, is presented in figure^J We can see that for those cases the region in the 
parameter space that corresponds to both conditions satisfied, (|12f) and (|13f) . and that 
gives the right relic density £lh 2 , is not so much extended with respect to the total 
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Figure 16. Experimental data (from )17| 'l confronted with PAMELA'S 
expectations for positron charge fraction for the DC model background. The 
propagation uncertainty band of the positron charge fraction and the curve that 
corresponds to the parameters of the best B/C fit in the middle are given for a 
better comparison. 



cosmologically allowed zone. The red zone mainly corresponds to small dumpiness 
factors, i.e. less than 20. 

In the tan (3 = 55 case (see the lower panel of figure I17|l there is the maximally 
extended cosmologically allowed region, and a huge part of it lies between the equi- 
clumpincss factors lines fd = 1 and fd = 10. So, the disentanglement of the 
supersymmetric signal is achieved without increasing the minimal dumpiness factors. 

In figure|3U]we gave two 3D plots that represent the functional dependance of the 
minimal fd factor from mo and m^: the left panel corresponds to figure IT8l while 
the right one corresponds to the lower panel of the figure El 

The last case is that of tan (3 = 60 (see figure 12 H for which the cosmologically 
allowed zone (red zone) is completely between the lines with fd=l and fd = 10. 

5.2. Some Examples of the Primary Component of the Antiproton Flux 

In this section we present different contributions to the antiproton flux induced by 
neutralino annihilations. Fluxes are calculated for different neutralino masses obtained 
from particular choices of the five mSUGRA parameters that satisfy the WMAP limits 
on the cosmological relic density and for different choices of the dumpiness factors fd. 

Higher neutralino masses produce primary contributions that improve high energy 
data fits but these models need higher dumpiness factors, due to the dependence from 
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Figure 17. Contour plots for the minimum fd needed for a PAMELA 
disentanglement (upper bounds of the translucent bands) and for the maximum 
fd allowed by current experimental data (lower bounds of the translucent bands). 
In the upper panel tan/3 = 50 while in the lower panel tan/3 = 55. The 
other parameters are Aq = and sgn(fi). Black color represents the regions 
in the parameter space that are excluded either by accelerator bounds or because 
electroweak symmetry breaking is not achieved or because the neutralino is not the 
lightest supersymmetric particle. The translucent regions denote the parameter 
space domains that correspond to models that satisfy both the conditions 1121 and 
fEfl . but now keeping fixed values of fd (those indicated inside the regions). Red 
(dark shaded) are domains with tth 2 in the WMAP region 0.09 < Qh 2 < 0.13, 
while green (light shaded) are the parameter space domains with 0.13 < Clh 2 < 
0.3. Under the line with fd = 1 there are no models that satisfy both the 
conditions 1121 and 1131 . We also show the equi- neutralino mass contours (blue 
dashed lines). 
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Figure 18. Left: equi-clumpiness factors lines as functions of mo and m 1 /2 
parameters in the first interesting part of the parameter space ("zoom" of the 
WMAP allowed zone) for tan = 50. Regions in the parameter space that are 
excluded either by accelerator bounds or because electroweak symmetry breaking 
is not achieved or because the neutralino is not the lightest supersymmetric 
particle are black. Domains with Qh 2 in the WMAP region 0.09 < Qh 2 < 0.13 
are red, while green indicates regions with 0.13 < Uh 2 < 0.3. Right: the same 
plot as in the left panel, but with blue region under the line with fd = 1 in which 
there are no models that satisfy the both conditions 1121 and 1131 simultaneously. 




Figure 19. Equi-clumpiness factors lines as functions of the mo and mi/2 
parameters in the second interesting part of the parameter space ("zoom" of 
the focus point region) for tan /3 = 50. Regions in the parameter space that are 
excluded either by accelerator bounds or because electroweak symmetry breaking 
is not achieved or because the neutralino is not the lightest supersymmetric 
particle are black. Domains with Qh 2 in the WMAP region 0.09 < Qh 2 < 0.13 
are red, while green indicates regions with 0.13 < Uh 2 < 0.3. 
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Figure 20. dumpiness factor as a function on the (mo, WI1/2) plane for fixed 
tan f3 = 50 in the "zoomed" zone from the figure [T%l Qeft panel) and for tan /3 = 55 
for all the masses we considered, i.e. the whole mass zone from the lower panel 
from the figure [T7I (right panel). 
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Figure 21. Same contour plot as in figure IT71 for tan j3 = 60, Ao = and 

sign(fi) = +1. 



the inverse neutralino mass squared m x 2 in the primary antiproton flux formula l|10|) . 

Small neutralino masses produce relatively high fluxes at small energies where 
the data are already fitted with a background contribution only. As a consequence a 
relatively small part of the mSUGRA parameter space gives the primary contributions 
that satisfies the experimental data fit. As can be seen in figure [22 there are models 
that give a total flux that is above the background uncertainty band for fd ~ 50 in 
the high energy range. We can consider two limiting classes of background models. In 
the first, most favorable class there are models near the lower limit of the uncertainty 
band. In this case it is always possible to add a supersymmetric contribution that 
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Figure 22. Different neutralino annihilations induced contributions to the 
total antiproton flux with the DC model background that corresponds to the 
propagation parameters of the best fit of the B/C data. The dashed lines 
correspond to the background contributions, the punctuated lines correspond 
to the neutralino induced contributions, while the solid lines correspond to 
total antiproton fluxes. As an example, in the lower panel we also show the 
PAMELA expectations for the total flux together with the uncertainty band for 
the background (translucent band between dot-dashed lines). All the three models 
are detectable by PAMELA. 
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gives fluxes compatible with the experimental data and detectable by PAMELA with 
a higher dumpiness factor. In the second, unfavourite class there are models near the 
upper bound of the uncertainty band. These models require much higher dumpiness 
factors and neutralino masses and this reduce the allowed zone in the supersymmetric 
parameter space. 

6. Conclusions 

In the first part of this work we deduced systematically the uncertainties of the most 
important cosmic rays spectra, with an emphasis on positrons and antiprotons. 

The DR and DRB models systematically overestimate the positron, protons and 
helium production. Electrons are overproduced at low energies, especially in the case 
with the break. The antiprotons are instead underproduced. For the model with 
diffusion and convection the results are in better agreement with the data. Further 
measurements of antiproton and positron spectra, primary to secondary CR ratios 
and solar modulation, as well as the precise determination and parametrization of 
important nuclear cross sections 03] seem to be crucial to determine the correct 
propagation model. 

In the DR and DRB models the antiproton experimental data fits can be easily 
improved adding different primary components coming from neutralino annihilations 
or from some other exotic contributions. In the framework of the DC model, exotic 
contributions remain possible at high energies (E > 20 GeV) and are not excluded at 
lower energies, due to the relatively large background uncertainties. 

In the second part of this work, we treated the detection of cosmic rays with the 
upcoming PAMELA experiment. PAMELA will measure with high statistics and in 
a wide energy range the spectra of protons, helium and light nuclei JSHI an d this will 
help to reduce the uncertainties in the antiproton and positron background fluxes. 
We computed the statistical errors for positrons and antiprotons (see figures figure IT2*1 
and figure IT3|) and we presented in figures ITK1 and IT^l PAMELA expectations for the 
antiproton proton ratio and positron charge fraction together with their respective 
propagation uncertainties. 

Then we studied the possibility to detect an eventual neutralino induced 
component in the antiproton spectra in the mSUGRA framework assuming the 
standard component to be the best fit for the DC model. In this case PAMELA 
will be able to disentangle a neutralino induced component for halo models that has 
fd as low as ~ 1. It is also possible to find models that give a total flux that is 
above the background uncertainty band but for larger value of fd (see lower panel of 
figure !^ . 

Background models near the lower limit of the uncertainty band still give fluxes 
compatible with the experimental data and detectable by PAMELA but for higher 
dumpiness factors. 

For background models near the upper bound much higher dumpiness factors 
and neutralino mass will be needed and this would reduce the allowed zone in the 
supersymmetric parameter space. 
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